EPJ manuscript No. 

(will be inserted by the editor) 



Dynamics of a polymer chain confined in a membrane 



o 

< 



o 

in 

e 

c 

o 
o 



S. Ramachandran 1 , S. Komura 1 a , K. Seki 2 , and G. Gompper 3 

1 Department of Chemistry, Graduate School of Science and Engineering, Tokyo Metropolitan University, Tokyo 192-0397, 
Japan 

2 National Institute of Advanced Industrial Science and Technology, Ibaraki 305-8565, Japan 

3 Institut fur Festkorperforschung, Forschungszentrum Jiilich, D-52425 Jiilich, Germany, EU 

Received: date / Revised version: date 

Abstract. We present a Brownian dynamics theory with full hydrodynamics (Stokesian dynamics) for a 
Gaussian polymer chain embedded in a liquid membrane which is surrounded by bulk solvent and walls. The 
mobility tensors are derived in Fourier space for the two geometries, namely, a free membrane embedded in 
a bulk fluid, and a membrane sandwiched by the two walls. Within the preaveraging approximation, a new 
expression for the diffusion coefficient of the polymer is obtained for the free membrane geometry. We also 
carry out a Rouse normal mode analysis to obtain the relaxation time and the dynamical structure factor. 
For large polymer size, both quantities show Zimm-like behavior in the free membrane case, whereas they 
are Rouse-like for the sandwiched membrane geometry. We use the scaling argument to discuss the effect 
of excluded volume interactions on the polymer relaxation time. 
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PACS. 82.35.Lr Physical properties of polymers 
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1 Introduction 

Integral membrane proteins play a vital role in a variety of 
cell functions such as solute transport, signal transduction 
and regulation of membrane composition [T] . Owing to fi- 
nite temperatures, such proteins along with other mem- 
brane components are constantly undergoing Brownian 
motions. The resulting diffusive motion plays an impor- 
tant role in determining their transport properties. Hence 
the studies of diffusion constitutes an important basis for 
understanding the physical properties of membrane pro- 
teins, and have been an active area of focus on model 
systems [2,3,4,5,6 as well as on living cells [71IB1IWTU] . 

Although proteins consist of polymeric units of amino 
acids, the standard approach is to consider them as rigid 
disks moving in a two-dimensional (2D) liquid membrane 
under low-Reynolds number conditions. The diffusion co- 
efficient of a rigid disk translating in a membrane which 
is embedded in a three-dimensional (3D) bulk fluid was 
calculated by Saffman and Delbriick (SD) [111112] . The 
obtained logarithmic size dependence is valid in the limit 
of small disk sizes. The SD theory was formally extended 
by Hughes et al. to all disk size ranges [12] . In the case of 
large disk sizes, they showed that the diffusion coefficient 
is inversely proportional to its size, which is analogous 
to the Stoke-Einstein relation in 3D. However, it should 
be noted that there is no single expression of the diffu- 
sion coefficient which covers the whole disk size ranges. In 
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a separate theoretical study, Evans and Sackmann (ES) 
employed a phenomenological approach to calculate the 
diffusion coefficient of a rigid disk moving in a membrane 
attached to a substrate [2] . The presence of a substrate in 
close proximity to the membrane was taken into account 
through a momentum decay term in the hydrodynamic 
equations. An extension of this work taking into account 
the effect of the advective terms has also been done |15j . 
In addition to these, diffusion of rod shaped objects on 
membranes |16pi7| or on Langmuir monolayers |18j have 
also been theoretically analyzed. 

Alternatively, the membrane protein can be regarded 
as a polymer chain rather than a rigid disk. In this case, 
the internal degrees of freedom of the polymer should be 
taken into account, which is the main subject of this paper. 
Apart from the protein analogy, hydrophobically modified 
polymers which adhere to the membrane could also be de- 
scribed using our description |19j . There are two theoret- 
ical works preceding the current work. One of them is a 
study by Muthukumar on the dynamics of a hydropho- 
bic polymer confined in a 2D liquid membrane |20) . In 
his theory, the membrane itself was treated as an isolated 
2D system having an anisotropic viscosity. It was shown 
that the mean squared displacement of a monomer obeys 
a diffusive law. He also pointed out that the mode de- 
pendence of the relaxation time arises from the excluded 
volume effect. The second and more direct precursor to the 
present work is the analytical calculation of the diffusion 
coefficient of a polymer chain |21) using a 2D hydrody- 
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Fig. 1. A transmembrane protein approximated as a poly- 
mer chain (red chain) embedded in a liquid membrane (yel- 
low particles). The membrane itself is surrounded by solvent 
(blue particles). Only a few representative solvent particles are 
shown. 



namic model with momentum decay [52,23,24,25 . Since 
their hydrodynamic model is essentially equivalent to that 
used by ES, the asymptotic size dependencies of the dif- 
fusion coefficient is the same between disks and polymers; 
namely, logarithmic in the small size limit, and algebraic 
in the large size limit. 

Further motivation is provided by the experiments on 
DNA molecules embedded on a cationic supported mem- 
brane [21)11271128] . The negative charge of the DNA molecules 
leads to strong adhesion with the membrane so that only 
the lateral motions are allowed. The measured diffusion 
coefficient showed a Rouse-like behavior. More recently, 
a similar experiment with DNA on a free standing mem- 
brane has been conducted [215]. Another related situation 
can be found in a dilute polymer solution confined be- 
tween narrow slits. Based on the scaling argument, such a 
polymer was predicted to show a Rouse-like behavior [301 
31,32,33 . In an attempt to verify these predictions, ex- 
periments on dilute solutions of DNA confined in narrow 
slits have shown that the exponents for conformation and 
chain relaxation of DNA is 2.2 which lies between 2D and 
3D behaviors [3~4] . 

In this paper, we discuss a Brownian dynamics theory 
for a polymer chain confined in a membrane. As schemat- 
ically presented in fig. [TJ we consider a polymer chain 
(connected red particles) embedded in a liquid membrane 
(yellow particles) surrounded by a 3D bulk fluid (blue par- 
ticles) and walls (not shown). We first derive the mobility 
tensors for the two geometries, i.e., a free membrane em- 
bedded in a bulk fluid, and a membrane sandwiched by 
two walls. For these two cases, we shall obtain the corre- 
sponding analytical expressions for the polymer diffusion 
coefficient which are valid for all sizes. We further perform 
a Rouse normal mode analysis to calculate the relaxation 
time and dynamical structure factor. For large polymer 
sizes, these quantities show Zimm-like behavior in the free 
membrane case, whereas they are Rouse-like in the pres- 
ence of walls. We also use a scaling theory to discuss the 
effect of excluded volume interactions on the relaxation 
times for the two geometries. The present work demon- 
strates the importance of the outer environment of the 



membrane in determining the dynamics of a 2D polymer 
chain. 

In the next section, we start to set up the governing 
equations for the membrane and the derivation of mo- 
bility tensors. With the introduction of the polymer, the 
general formalism of the problem is constructed in sect. [3] 
Sections @] and [5] describe the results for the polymer dy- 
namics for the free and confined membrane limiting cases, 
respectively. Excluded volume effects are discussed using 
scaling arguments in sect. [SI We finally close with several 
discussions in sect. [7] 



2 Membrane hydrodynamics 

Before introducing the polymer, we first establish the gov- 
erning equations for the membrane and its surrounding 
environment. The aim of this section is to derive the mem- 
brane mobility tensors which will be used in the later sec- 
tions for the polymer equations of motion. The present 
calculation closely follows the formulation by Inaura and 
Fujitani [35) . However we consider a more general situ- 
ation which will be described below. The details of the 
calculation are relegated to appendix A. 

As shown in fig. [2J we assume that the membrane is 
an infinite planar sheet of liquid, and its out-of-plane fluc- 
tuations are totally neglected, which is justified for typi- 
cal bending rigidities of bilayers. Relaxation dynamics of 
membrane fluctuations near walls have been previously 
considered [36,37,38 . The liquid membrane is embedded 
in a bulk fluid such as water or solvent which is bounded 
by hard walls. Let v(r) be the 2D velocity of the mem- 
brane fluid and the 2D vector r — (x,y) represents a point 
in the plane of the membrane. We first assume the mem- 
brane to be incompressible 

V-v = 0, (1) 

where V is a 2D differential operator. We work in the low- 
Reynolds number regime of the membrane hydrodynamics 
so that the inertial effects can be neglected. This allows 
us to use the 2D Stokes equation given by 

r/V 2 v - Vp + f s + F = 0, (2) 

where r\ is the 2D membrane viscosity, p(r) the 2D in- 
plane pressure, f s (r) the force exerted on the membrane 
by the surrounding fluid ("s" stands for the solvent), and 
F(r) is any other force acting on the membrane such as 
that due to a polymer chain introduced in sect. [3[ 

As presented in fig. [2] the membrane is fixed in the xy- 
plane at z = 0. The upper (z > 0) and the lower (z < 0) 
fluid regions are denoted by " + " and " — " , respectively. 
The velocities and pressures in these regions are written as 
v ± (r,z) andp ± (r, z), respectively. Since the 3D viscosity 
of the upper and the lower solvent can be different, we 
denote them as 77^ respectively. Consider the situation in 
which impenetrable walls are located at z — ±h , where 
h + and h~ can be different in general. We note that this 
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Fig. 2. Schematic picture showing a planar liquid membrane 
having 2D viscosity r\ located at z — 0. ft is sandwiched by a 
solvent of 3D viscosity 77* . Two impenetrable walls are located 
at z — ±/i bounding the solvent. 



point differs from ref. |35j . Similar to the liquid membrane, 
the solvent in both regions are taken to be incompressible 



0, 



(3) 



where V represents a 3D differential operator. We also 
neglect the solvent inertia and hence the solvent obeys 
the 3D Stokes equations 



?7 s ± V 2 v ± - Vp ± = 0. 



(4) 



The presence of the surrounding solvent is important be- 
cause it exerts force on the liquid membrane. This force, 
indicated as f s in eq. ([5]), is given by the projection of 
(er+ — <t~) z= o ■ e z on the xy-plane. Here e z is the unit 
vector along the z-axis, and a ± are the stress tensors due 
to the solvent 



<J = -J^I + T^Vv* + (Vv*) 1 ] 



(5) 



In the above, I is the identity tensor and the superscript 
"T" indicates the transpose. 

Using the stick boundary conditions at z — and z = 
zt^, we solve the hydrodynamic equations © and ^ to 
obtain f s . Then we calculate the membrane velocity from 
eq. @ as 

v[k]=G[k]-F[k], (6) 

where v[k] and F[k] are the Fourier components of v(r) 
and F(r) defined by 



and 



v(r) = 



F(r) 



/(l^ V[k]eXp( * k ' r) ' 



d 2 k 

(27T) 2 



F[k] exp(ik • r), 



(7) 



(8) 



respectively, and k = {k x , k y ). After some calculations (see 
appendix A for the details), one can show that the mobil- 
ity tensor G[k] in Fourier space is given by 



G Q0 [k] 



1 



r\k 2 



k[r/t coth(kh + ) 

k a kji 



T] s coth(fc/i~)] 



SaB — 



k 2 



(9) 



with a, (3 = x,y and k = |k|. As in ref. [35], we mainly 
consider (except in sect. [7]) the case when the two walls 
are located at equal distances from the membrane, i.e., 
h + = h~ = h. Then the above mobility tensor becomes 



G a/3 [k] = 



1 



77 [k 2 



vk coth(kh)] 



s a 







k a k 



k 2 



(10) 



r) s )/2. An almost equiv- 



where v = 2r) s /r] with r) s 
alent expression to eq. (fTU|) has also been derived for Lang- 
muir monolayers in which there is only one wall or a sub- 
strate [18, 39 . In the following, we discuss the two limiting 
situations of eq. (fTt)|) . 

Saffman and Dclbruck (SD) investigated the case when 
the two walls are located infinitely away from the mem- 
brane pTTll 1 2 j . This is called as the free membrane case 
and all the related physical quantities are denoted by the 
superscript "SD". Taking the limit of kh ^> 1 in eq. (fTU|) . 
the mobility tensor becomes [35l [40ll4Tj 



Ok] = 



1 



rj{k 2 + vk) 



5at 



k a kj^ 



(11) 



Notice that the length v~ x is called the SD hydrodynamic 
screening length. The real space expression of this mobility 
tensor is obtained by the Fourier transform of eq. (jlll) 



G(r) 



d 2 k 
(2n) 2 



G[k] exp(ik-r). 



(12) 



Performing the calculations presented in appendix B, we 
obtain I51I40U41I 



G^(r) 



1 

477 



H (i^)-Fo(w 
HiH , Yi (vr) ' 







1 

4?7 / KV 2 r 2 
2Y 1 (vr) 



vr 
vr 

H (vr) + Y a (vr) 



(13) 



where r = |r|. In the above, H„(z) are the Struve func- 
tions and Y n (z) are the Neumann functions or the Bessel 
functions of the second kind. 

In the opposite kh <C 1 limit, the membrane is con- 
fined between the two walls. Since such a limiting case 
was considered by Evans and Sackmann (ES) [11], we de- 
note all the physical quantities for this situation with the 
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superscript "ES" . In this case, eq. (jlOp takes the following 
form 



u a/3 



[k] = 



1 



Sal: 



k a k 



k 2 



(14) 



T](k 2 + K 2 ) 

where k = (y /h) 1 / 2 . This new length scale k _1 is the ES 
hydrodynamic screening length. We note that k _1 is the 
geometric mean of v~ x and h [32]. The above ES mo- 
bility tensor was previously used in a phenomenological 
membrane hydrodynamic model 21,24,25,41,43 . Follow- 
ing the calculations in appendix B, the real space repre- 
sentation of the ES mobility tensor becomes 



1 

2lXT\ 
1 

27T77 



K (nr) -+ 
-K (K,r) 



2K 1 ( K r) 



1 



?a/3 



r ar/3 



(15) 



where K n (z) are the modified Bessel functions of the sec- 
ond kind. We use either eq. (fTTj) or eq. ([T3]l in our subse- 
quent calculations with a polymer. 

Strictly speaking, eqs. (fK3"]l and ([T5"]) should have been 
obtained by taking the limits of vh 3> 1 and <c 1, 
respectively, after the integration of eq. (fit)]) over k. How- 
ever the inverse Fourier transform of eq. (|10p is nontrivial. 
The present approach serves our purpose and the rigorous 
derivation will be given in a separate publication. 



3 Dynamics of a 2D Gaussian polymer chain 
embedded in a membrane 

We are now ready to introduce a polymer into the mem- 
brane. For simplicity, we first work with a Gaussian poly- 
mer chain whose conformation is given by a set of N posi- 
tion vectors denoted as {R„} = (Ri, . . . , Rat) embedded 
in the 2D membrane. The excluded volume effects will be 
discussed later in sect. [6] It is implicitly assumed that the 
polymer relaxation time is much longer than that for the 
typical hydrodynamic disturbances so that the membrane 
can be effectively considered as a 2D liquid. If the polymer 
consists of monomers which exert a set of point forces f„ 
acting at R„, the external force due to the polymer F(r) 
in eq. ([2]) can be written as 



P(r) = ^f„5(r-R n ). 



(16) 



In writing this expression, we have assumed that the su- 
perposition principle holds. With the use of the mobility 
tensor obtained in the previous section, eq. (J2J can be 
formally solved as 



v(r)=^G(r-R n ).f„ 



(17) 



Since monomers move with the same velocity as the mem- 
brane, their velocities are given by 



dRn(t) 

dt 



— v (R?i) — ^ ' G JlTO • f„ 



(18) 



where we have used the notation G nm = G(R„ — R m ). 

The Langevin equation for a polymer chain embedded 
in a membrane is written as [34,45 



dt 



dU 

nm ' \ ~ ;, D T SmrJ 
Oti. m 



kBT^_d_ 

2 2^dR, „, " 



(19) 



where Cm{t) is the Gaussian random force acting at R m , 
k B the Boltzmann constant, T the temperature. The po- 
tential energy of the 2D polymer has the form 



U 



k B T 
b 2 



^~^(Rn — Rn-l) 2 , 



(20) 



n=2 



where b is the Kuhn length. It can be shown that the 
mobility tensors eqs. (1331 and (IT5|) satisfy 



d 



<9R„ 



• G nm — 0. 



(21) 



Substituting eq. ([20]) into eq. ([T9"]) and using eq. (|2Tj) . we 

have 



dR n (t) 
dt 



5> 



2k B Td 2 R m (t) 
b 2 dm 2 



+ Cm(t)). (22) 



Due to the hydrodynamic coupling between different 
parts of the polymer, the above equation is non-linear and 
difficult to solve analytically. In order to overcome this dif- 
ficulty, we employ the preaveraging approximation which 
has been successfully used for a polymer in 3D solvent [451 
146] . Assuming that the polymer is close to its equilibrium, 
we replace G nm by its equilibrium value (G nm ) such that 



(G r , 



d{R„}G„ m *({R„}) 



1 



■ cxp 



/,, ir\n — m\b 2 ^ \ \n — m\b 2 

= g(n - m)I, 



G nm (r) 
(23) 



where $({R„}) is the 2D Gaussian distribution function. 
Within this approximation, eq. (|22|) can be simplified as 

dK n (t) ^ f 2k B T d 2 R m (t) \ 

-- 2^9{n-m) ( — KZ^—+Cm(t) ) . (24) 



Dt 



\ b 2 dm 2 



The above equation can be rewritten in terms of the 
Rouse normal coordinates defined by [35] 

i r N 

XpW " I d " cos (^r) R «W> (25) 



as 



^ = Eft.[-W)+c.w] 



(26) 
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with 



4TT 2 k B T 

Nb 2 



-P ■ 



(27) 



Here p, q = 0,1,2,... and g pq is the mobility tensor in 
terms of the normal coordinates: 



fqTrm\ 

\-^)9{n-m). 



f N dn f dm fpim\ 
9 ™ = J iv/o -iV C ° S (— J cos 

(28) 

If one neglects the contribution from the off-diagonal com- 
ponents of g pq , we finally obtain 



axp(t) 
dt 



g p [-kpKp(t) + C P (t)} , 



(29) 



with the definition g p = g pp . The Gaussian random forces 
Cp(t) satisfy the following conditions 

(C pa (t))=0, (30) 
{C pa {t)C q p{t')) = 28 pq 5 a p{g p )- 1 k B T5{t - t'). (31) 

Therefore the relaxation time of a polymer in terms of the 
Rouse modes is given by 



(32) 



g p k p 



Furthermore, the polymer diffusion coefficient can be cal- 
culated according to the following equation 



D = k B Tg = k B T^ ^ ^9(n-m) 



(33) 



Another useful quantity that we calculate is the dy- 
namic structure factor defined by 



(34) 



Since R„(i) — R m (0) is a linear function of Cn(t) obey- 
ing the Gaussian distribution, the distribution of R„ (t) — 
R m (0) is also Gaussian [47]. Hence we have 



(exp(ik- [R„(0-R™(0)])) 

= exp^-^((R„(t)-R m (0)) 2 ; 



(35) 



in 2D. Denoting the center of mass by Xo and using the 
inverse relation of eq. (|2"5)) 



R„ — Xo + 2 X 



P cos(— j 



(36) 



we can calculate the dynamical structure factor as 



S(M) 



l 

N 



J] exp 



k 2 Dt- ^\n-m\b 2 k 2 



Nb 2 k 2 ^ 1 fpirn\ /pirm 

p=l 1 



N 



[1 - exp(-t/T p )] 
(37) 



This completes the general formalism of the 2D polymer 
dynamics confined in a liquid membrane. In the next sec- 
tion, we consider the free and confined membrane cases 
separately by using eqs. (fTl"]) or ([T4"]) for the mobility ten- 



4 Polymer dynamics: free membrane case 

When the two walls in fig.[5]are located at infinite distance 
from the membrane, one can use the SD mobility tensor 
given by eq. (TlTj) . We calculate the preaveraged mobility 
tensor, relaxation time, diffusion coefficient and structure 
factor following the recipe described in the previous sec- 
tion. 



4.1 Mobility tensor 

In order to perform the preaveraging of the mobility ten- 
sor, we take the configurational average of eq. (jlll) by 
using the equilibrium probability distribution function of 
a Gaussian polymer in 2D: 



(G 



SD < 



d 2 k I - kk 

(2tt) 2 T](k 2 + vk) 



exp[ik • (R„ - R m )]y , 
(38) 



where k denotes a unit vector along k. This leads to 
(G™)=g s °(n-m)I with 



g SD (n-m) 

1 f°° k (I 

= - — / dk —z — exp f -b 2 k 2 \fi m 

47777 J k 2 + vk \ 4 1 

= exp ( — b 2 v 2 \n — ml 

87T/7 P V 4 1 

1, 



7rerfi ( —bi/y/\n — m\ 



Ei | -b 2 v 2 \n-m\ 



where erfi(z) is the imaginary error function 

erfi(z) = — ierf(iz), 

with erf (z) being the error function 

2 f z 2 
erf (2) = — = / due~ u , 



whereas Ei(— z) is the exponential integral given by 



Ei(-z) 



du ■ 



(39) 



(40) 



(41) 



(42) 



It should be noted that the obtained g sr> (n — m) is real 
despite the presence of complex functions. 
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4.2 Relaxation time 



In order to obtain the polymer relaxation time, we first 
substitute eq. ([3"51 into eq. (f^5|) to express the mobility 
tensor in terms of the Rouse normal coordinates 



„SD 



1 



TTTjNb 2 

n z p - 



dk 

(k 2 

^l 2 5 - 



fc 3 



vk)[k 4 + (4irp/Nb 2 y 



21n(7rp/ S 2 )S 2 + (VWp)* 3 



16irrj 



ir V + S 4 



(43) 



(note that <7p D = g^p). In the above, we have defined 

the dimensionless polymer size 8 = \J~Nbvj2. Since the 
radius of gyration for the 2D Gaussian polymer is R g = 
\J~Nbj2 [45], S can be also written as <5 — R g v. Using 
eq. (I3"2"j) . the relaxation time becomes 



SD 



ANb 2 r) 
■nk B T 



7T 2 P 2 



p 2 [ir 2 p - y/2p~n 3 / 2 5 + 2 ln(TT P /S 2 )S 2 + (VWp)<5 3 ] 



(44) 



This expression shows how the presence of the bulk solvent 
affects the relaxation time. 

We consider two asymptotic limits of eq. (14"4"1) . For 
small polymer sizes or 5 < 1, we have 



-SD 



ANb 2 i] 1 
itk^T p 



(45) 



For large sizes, the condition 8^1 yields 



SD 



ANb 2 f] 



irk B T V27rp 3 / 2 



(46) 



It should be noticed that eq. (|4*S1) depends only on the 
solvent viscosity r] s but not on the membrane viscosity r\. 

In fig. [21 the scaled relaxation time eq. (|44l) is plotted 
as a function of S for p — 1 and 10 as solid lines. For 
small (5, the relaxation time is independent of the polymer 
size, which is consistent with eq. f)45|) . The relaxation time 
increases through a crossover regime towards a linear be- 
havior as given by eq. (|4"oT) . Such a crossover occurs around 
the region where the polymer size R g is comparable to the 
SD hydrodynamic screening length v~ x , i.e., 6 ~ 1. In the 
limit of large S, the p-dependence of the relaxation time 
is analogous to that obtained from the Zimm model [4"5] . 
The solid lines in fig. 2] show the relaxation time as a 
function of the Rouse normal mode p for S — 0.01 and 
100. These lines have slopes —1 and —3/2 for 6 — 0.01 
and 100, respectively. These results indicate the different 
mode dependencies in the two limiting polymer sizes. 




Fig. 3. Scaled relaxation time irksTrp / 'ANb 2 "r\ as a function 
of 5 = RgV for free membranes (solid lines) or e = R g K for 
confined membranes (dashed lines) for p = 1 and 10. 




Fig. 4. Scaled relaxation time irk^TT p / ANb 2 rj as a function 
of p for different values of 5 = 0.01, 100 (solid lines) and 
e — 0.01, 100 (dashed lines). The two curves for 5 = e = 0.01 
overlap each other and cannot be distinguished. 



4.3 Diffusion coefficient 

By substituting eq. (I3TR) into eq. ([3"3")l . the diffusion coeffi- 
cient of the polymer can be obtained as 



D 



SD 



k B T 1 



(7rerfi(<5) - Ei(<5 2 )) exp(-,5 2 ) + 



4v^ 



47r?7 S 4 

+ 5 2 - {\nS 2 + j)(5 2 - 1) - 2^8 



(47) 



where 7 = 0.5772 • • • is Euler's constant. This expression 
for the polymer diffusion coefficient is valid for all the 
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Fig. 5. Dimensionless diffusion coefficient Anr/D/kBT as a 
function of S = R s v (solid line) or e = R g n (dashed line). 



ranges of <5. Equation (l47l is one of the main results of 
this paper. 

We now discuss the asymptotic limits of eq. (14"T1) for 
small and large polymer sizes. When 8 <C 1, it reduces to 



D 



SD 



47T77 V 24 



Such a logarithmic behavior is consistent with that of an 
object in a pure 2D system [TT1H2] . In the opposite limit 
of 8 ^> 1, we have 



D 



SD 



Atti] 38 6y/irr/ s R g 



Similar to eq. (|46[) . this expression depends only on r) s . The 
obtained l/?7 s i?g-dependence is analogous to that of an 
object moving in 3D fluid as well as the result by Hughes 
et al. [T3]. 

In fig. [5j we plot the diffusion coefficient D as a 
function of 8 (solid curve). With the increase in the poly- 
mer size, there is a crossover from logarithmic to algebraic 
decay indicated by eqs. (14"51) and (14"9"1) . respectively. The 
change in the behavior of a domain of the diffusion coef- 
ficient with the addition of solvent has also been shown 
through recent dissipative particle dynamics simulations 
on pure 2D and quasi-2D systems [15] . 



4.4 Dynamic structure factor 

The last quantity calculated for the free membrane geom- 
etry is the dynamical structure factor S SB (k, t) defined in 
eq. ([3"4"]1 . Since the full expression is rather complicated, 
we derive several analytical expressions for the limiting 
cases. For R g k <C 1, we have 



S bU (k,t) w Nexp(-k 2 D bU t) 



(50) 



where D is given by eq. (jT7]) . This is reasonable because 
only the center of mass motion of the polymer is captured 
in the small angle regime. 

For R g k ^> 1, on the other hand, we only consider the 
time region t <C r| D because S SB (k, t) becomes very small 
for t ^> r^ D . In this case, we have 



S SD (M)«l]Texp 



1, 



2;„2 



m\b k 



Nb 2 k 2 



E 

P =i 



L fpirn\ 

2 cos VW ) cos 



b 2 k 2 



du exp ( — u — k 2 I\{u)), 



(51) 



with 
h(u) 



Nb 2 
2^" 



dp — cos 



( 4irpu 
\Nb 2 k 2 



-exp(-Vr p SD )]. 

(52) 

When <5 3> 1, we can use the limiting expression for t~ d 
as obtained in eq. ([4"6"]) . In this case, the above expression 
becomes 

S su (k,t) 
8 



b 2 k 2 



du exp 



-u~(r SD t) 2 / 3 Wl ((r SD ty 2 ^u) 



(48) with the decay rate 



and 



(49) w i( u ) 



dx 



k B Ty/Nbk 3 _ k B Tk 3 
l6i]S 8r]u 

cos(xu) 



(53) 



(54) 



exp 



(-x 3 / 2 /V2)]. (55) 



Notice the fc 3 -dependence of the decay rate. For r° B t 3> 1, 
the above expression is further simplified to 

S SD (k, t) a 5 SD (k, 0) exp(-1.35(r SD i) 2 / 3 ), (56) 

since u>i(0) = _T(l/3) ~ 1.35. This expression is valid in 
the limit of large polymer sizes such that R g k 3> 1 and 
<5> 1. Large wave vectors (probing the internal motion of 
the polymer) and hydrodynamic screening lengths (high 
solvent viscosity) will also lead to the same expression. 
The applicable time window for eq. ([55]) is 1/ r su <C t <C 
Tp D . These expressions for large 8 are analogous to that 
obtained from the Zimm model [45 in the R g k 3> 1 regime 
for polymer in 3D. 



5 Polymer dynamics: confined membrane case 

When the thickness of the solvent layers is very small, the 
membrane is now almost confined by the two walls. How- 
ever, there is a thin lubricating layer between the mem- 
brane and the walls so that h ^ 0. In this case, we use the 
ES mobility tensor given by eq. (|14[) . Several quantities 
for this case are obtained below. 
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5.1 Mobility tensor 



5.3 Diffusion coefficient 



By using eq. (O, the preaveraged mobility tensor is cal- 
culated from 



/ f d 2 k I — kk 

(G™) = I j (^I^^y ex P [ lk • (R„ - R m) ] 



(57) 



This results in (G E ,S) = g ES (n ~ m)I with 

g ES (n - to) 

1 f°° k ( 1 

= / dk — - exp ( — b 2 k 2 \n — to 

4tt77 J k 2 + k 2 1 \ 4 1 

= exp ( — b 2 K 2 \n — ml] Ei ( — b 2 K 2 \n — m 

8nri V 4 7 V 4 



which was previously derived in ref. |21j 



(58) 



5.2 Relaxation time 

The polymer relaxation time can be obtained by substi- 
tuting eq. ([551) into eq. (T2"gj) . Then we have 



i 



dfc- 



1 7r 2 p + 2e 2 ln(£ 2 /(7rp)) 

167T77 



(59) 



In the above, we have defined the dimensionless polymer 
size as e = \fNbn/2 = R g n which should be distinguished 
from S in the previous section. Then the relaxation time 
can be written as 



ES 



4Nb 2 7i 



n 2 p 2 



irk B T p 2 [ir 2 p + 2e 2 ln(e 2 /(7rp))] ' 



In the limit of e <C 1, it reduces to 

t es w 4Nb 2 v 1 
p nk^T p' 



(60) 



(61) 



which coincides with eq. (|45|) . In the opposite limit of e 
1, one gets 



4Nb 2 T] 



irk B T 2p 2 ln(e 2 /(7rp))' 



(62) 



In fig. [3l we plot t es as a function of e for p — 1 and 10 
in dashed lines. The algebraic e 2 -dependence in eq. (l62|) is 
seen for large e. The dashed lines in fig. 0] are the plots of 
Tp S as a function of p for e = 0.01 and 100. For e = 100, 
the slope —2 is consistent with eq. (j52"j) neglecting the 
logarithmic correction. Notice that this p-dependence is 
in contrast to that for the free membrane case given in 
eq. gSJ. 



With the use of eq. (|58|> . the diffusion coefficient for the 
confined membrane geometry is written as 



D ES = 



fc R T 1 



[(l+e 2 )(21n £ + 7 ) 



47r?7 e 4 
- exp(e 2 )Ei(-e 2 )]. 

This equation was also obtained before 
expression for e <C 1 is 



D 



ES 



k B T 



lne - 



7 



(63) 

The limiting 



(64) 



which coincides with eq. (|48|) as long as e is replaced by 
5. When e> 1, eq. (|53l reduces to 



ES 



fc B T 1 



k B Th 



Attt] e 2 8irr] s Rl 



(65) 



This l/i? 2 -dependence is a characteristic of a system in 
which there is momentum loss from the membrane to the 
surrounding environment |50) . An intuitive understanding 
of the large size behaviors of the diffusion coefficient in 
terms of the conservation principles will be described in 
sect. [3 The dashed line in fig. [5] shows the plot of D ES 
as a function of e, showing the logarithmic and algebraic 
behaviors as derived above. 

According to the definition of 5 and e, one obtains 
e = S/y/uh. This correspondence leads to a rescaling of 
the e-axis in figs. [3] and [5] Here the value of y/vh cannot 
be taken arbitrarily since the condition for the confined 
membrane is given by vh <C 1 as explained in sect. [5] 



5.4 Dynamic structure factor 

The dynamic structure factor can be calculated in the 
same manner as before. For R g k <C 1, we have 



S^foi) w A^exp(-fc 2 J D ii ^) 
For R g k > 1 and t < r ES , we get 



(66) 



S ES (k,t) 



b 2 k 2 J 



du exp( — ii — k 2 I 2 (u)), (67) 



with 



= "oZT / dp-, cos [1 - exp(-t/rf)]. 



2tt 2 



p 2 \Nb 2 k< 



(68) 



Considering £> 1 and neglecting the logarithmic depen- 
dence in eq. (|62[) . the above expression becomes 



b 2 k 2 



du exp -u- (r ES t)^ 2 w 2 ((r^t)-^ 2 u) 



(69) 
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with the decay rate 

_, ES _ k B TNb 2 k 4 _ k B Tk A 
32irrje 2 8nr]K 2 ' 

cos(xu) 



and 



w 2 (u) 



dx 



-[l-cxp(-x 2 )]. 



(70) 



(71) 



Note that r ES is proportional to fc 4 as for the Rouse 
model. For T ES t > 1, it can be further simplified to 

S ES (k, t) « S ES (k, 0) exp(-1.13( J T ES i) 1/2 ), (72) 

since ^2(0) = 2/y/n f=s 1.13. This expression is valid when 
R g k 1 and e > 1. The relevant time interval for this 
expression is l/T ES < t < r ES . 



6 Excluded volume effects 

So far, we have treated only a 2D Gaussian polymer chain. 
In this section, we briefly discuss the effects of excluded 
volume on the dynamical quantities. Even if the polymer 
is not a Gaussian chain, one can still use the preaverag- 
ing approximation for the mobility tensor. For the free 
membrane case, one can generally show that 



1 



g (n-m) = —(H (vr nm ) - Y (ur nm )), 
877 



(73) 



where r nrn = |R„ — R m | (see appendix B). Then the dif- 
fusion coefficient can be expressed as 



D 



SD _ B 



k B T 



81, 



f N dn f N dm 

J 77 J "]y"( Ho( ^ rnm ) - Y o(vr nm )}. 



Similarly for the confined membrane case, we find 



D 



ES 



k B T 
4irn 



f N dn f N dm 
Jo nJ q lf( K ^ m )), 



(74) 



(75) 



(see eq. (|B.14[) ). These expressions are the 2D analog of 
the Kirkwood formula 45 . It should be emphasized that 
they are rigorous even in the presence of excluded volume 
effect. 

For excluded volume chains, however, the appropriate 
equilibrium distribution function $({R„}) needed to cal- 
culate averages is not known [55]. Therefore, we cannot 
obtain rigorous forms of diffusion coefficient. Instead, we 
shall make use of scaling arguments to infer the effects of 
excluded volume interactions. Here, we limit our discus- 
sion to small and large polymer sizes. A simple argument 
is that the excluded volume effects lead to a rescaled poly- 
mer radius of gyration R g . Within the Flory theory, the 
radius of gyration scales as R g ~ bN Vw using the Flory 
exponent up. Up to a numerical factor, the diffusion co- 
efficient in the limiting cases will still show the same size 
dependence as in eqs. (14*51) . (|4*^1) . ([Ml) and (f6"5"]) in which R g 
is now replaced with that of excluded volume chains. In 



Table 1. Comparison of the relaxation times between a Gaus- 
sian polymer chain and a polymer with excluded volume inter- 
actions. 



limits 


ideal chain (up = 1/2) 


real chain (uf = 3/4) 


<5< 1 




T SD_ p -3/2 


£ < 1 




T ES ~p" 3 / 2 


8 > 1 






£ > 1 


r ES ~p- 2 


r ES ~ p" 3 



other words, we replace N with N 2vf in these equations. 
On the other hand, the relaxation time in the presence of 
the excluded volume effects can be obtained by replacing 
N/p with (N /p) 2vF in the Gaussian chain cases. In the 
small size limits, i.e, S, e <C 1, we have 



and 



SD 



_ES 



■qb 2 
kWT 



N\ 2V¥ 
P J 



r]b 2 fN 



k B T \ p 



2// F 



(76) 



(77) 



In the opposite of large polymer size limits, i.e., 5, e 3> 1 
we get 

k B T\ P ) ' 



SD 



(78) 



and 



T]b 4 K 2 

k B T 



P J 



(79) 



These are the scaling predictions. Notice that the above 
results can also be obtained by using the relation r ~ 
R 2 /D where R g includes the excluded volume effect [ST] , 
Table Q] shows the comparison between a Gaussian 
chain polymer (up = 1/2) and a chain with excluded vol- 
ume interactions (up = 3/4) in 2D. The former Gaussian 
case recovers all the relations in the previous sections (see 
eqs. @5J), (EED, @BJ), {52])). For the chain with up = 3/4, 
the relaxation time shows a p -3 / 2 -dependence for both 
the free and confined membrane cases when 5, e <C 1 . This 
is consistent with the result in ref. [3D]. For i > 1 and 
e ^> 1, we have r SD ~ p~ 9 / 4 and r ES ~ p~ 3 , respectively. 
These exponents are unique for polymers confined in a 
membrane. 



7 Discussion 

In this paper, we have investigated the dynamics of a 
Gaussian polymer chain confined in a liquid membrane 
taking into account the surrounding environment. For the 
most general geometry with the membrane, solvent and 
walls (see fig. [2]), we have derived the mobility tensor in 
eq. ([9]). We obtained the analytical expressions for the re- 
laxation time, diffusion coefficient and dynamic structure 
factor of a Gaussian chain for the two limiting cases of 
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free and confined membranes. These quantities were cal- 
culated within the preaveraging approximation of the mo- 
bility tensor in order to avoid the non-linearity introduced 
by the hydrodynamic coupling. We shall summarize and 
discuss the results obtained in this paper. 

Our theoretical analysis relies on the preaveraging ap- 
proximation, which is used to decouple the polymer fluc- 
tuations from hydrodynamics. The validity of this approx- 
imation has been studied in some detail for polymers in 
a 3D bulk fluid, where the preaveraging approximation 
(Zimm model) yields results which are not very different 
from more sophisticated calculations [55]. The response 
of a dilute polymer solution to an external field has also 
been experimentally verified to follow the predictions of 
the Zimm model |52) . Deviations from the Zimm model 
have been attributed to non-Gaussian distributions, the 
slowness to reach asymptotic behaviors or the effect of hy- 
drodynamic fluctuations, although a clear conclusion has 
not been reached |53j . Up to now, there have not been any 
experimental studies of the dynamics of polymers confined 
to membranes, to the best of our knowledge. Experiments 
using a combination of Langmuir trough as well as light 
scattering techniques should be able to test our predic- 
tions. The success of preaveraging in three-dimensional 
polymer solutions encourages us to expect the same to 
hold for quasi-2D systems. 

We first compare the relaxation times between the two 
cases. As seen from eqs. (l45|) and (joT]) for small polymer 
sizes, i.e., 6, e <C 1, both Tr D and t e show exactly the 
same mode dependence. This is due to the fact that for 
polymer sizes smaller than than the hydrodynamic screen- 
ing lengths (v^ 1 or the outer environment does not 
affect the polymer dynamics. The behavior r^ ~ p~ 3 / 2 
for 5 3> 1 (see eq. is analogous to the Zimm relax- 
ation time of a Gaussian polymer in a 3D solvent. On the 
other hand, the dependence t e ~ p~ 2 for e >• 1 (see 
eq. (fl)2")) ) is similar to the Rouse relaxation time in 3D. 

One of the important results of this paper is the deriva- 
tion of eq. (|4"7)l which provides the diffusion coefficient 
valid for all size ranges in the free membrane case. Accord- 
ing to eqs. igSJ) and (fS4"]) . the diffusion coefficients D SB 
and D ES show the same logarithmic dependence when 
S, e <C 1. Again this occurs when the polymer size is 
smaller than v~ l or A crossover from logarithmic to 
algebraic behavior takes place when the polymer size R g 
is comparable to these hydrodynamic screening lengths. 
Thus v^ 1 or kT 1 determines the length scale at which the 
outer environment surrounding the membrane becomes 
important. For a pure 2D system, the screening length 
is infinite, which causes the Stokes paradox. 

When the polymer size becomes much larger than the 
hydrodynamic screening lengths, the interactions are no 
longer only through the membrane. In the free membrane 
case, the outer fluid plays the role of the mediator of the 
hydrodynamic interactions. This attributes a 3D nature 
to the polymer dynamics and hence the scaling Z) SD ~ 
l/r/ s i? g is recovered (see eq. (|49]) ). In the confined mem- 
brane case, the presence of the walls takes away the mo- 
mentum from the membrane. Owing to the stick bound- 



ary condition, a linear shear velocity gradient is set up in 
the intervening fluid layer between the membrane and the 
walls leading to a momentum leakage from the membrane. 
As given in eq. ([55)1 , the resultant diffusion coefficient be- 
haves as D ES ~ l/i2g. These asymptotic behaviors are 
consistent with the results of correlated diffusion obtained 
in ref . [41] . 

The different large size behaviors of the diffusion coef- 
ficient can be better understood by focusing on the mobil- 
ity tensors. It is essentially described as a consequence of 
the conservation of mass and momentum principles |50J. 
The mobility tensor G(r) acts as a Green's function re- 
lating a source of disturbance at the origin to the velocity 
at a point r. At sufficiently large distances, the source 
of disturbance in the fluid can be thought of as a force 
monopole which introduces a source for momentum in the 
fluid. A moving object also causes a perturbation in the 
mass density and this can be regarded as a mass dipole 
(a source and a sink of mass density). We now apply the 
conservations of momentum and mass to the 3D and 2D 
cases separately. 

For the 3D case, the presence of a force monopole at 
the origin should cause the momentum flux or stress a de- 
cay as 1/r 2 in order to conserve the total momentum (r 2 
being proportional to the area of a sphere surrounding 
the force monopole). Since the shear stress is related to 
the fluid velocity through a ~ rj s v/r, we have v ~ l/r/ s r. 
This implies that the mobility tensor should also scale as 
l/r/ s r. Concerning the mass conservation principle, a mass 
monopole would create a flow velocity that decays as 1 /r 2 . 
Since we have a mass dipole, the resulting velocity now de- 
cays as 1/r 3 . Comparing these two effects, the contribu- 
tion to the velocity due to momentum conservation (which 
varies as 1/r) always dominates at large distances in 3D. 
This essentially explains the behavior D su ~ l/rj s R g . 

Let us now consider the 2D case in which the mem- 
brane is in contact with the walls leading to a loss of 
momentum from the membrane. This implies that the 
momentum is not conserved, and the only contribution 
to the velocity is from the mass conservation. In 2D, a 
mass monopole will create a velocity which decays as 1 /r 
(r being the perimeter of a circle surrounding the mass 
monopole). Hence the velocity and the mobility tensor due 
to a mass dipole decays as 1 /r 2 . This explains the scaling 
D ES ~ 1/i? 2 . These two different behaviors of the mo- 
bility tensors are reflected in the diffusion coefficients for 
large size polymers. Incidentally, for the pure 2D case, the 
stress decays as 1/r due to the momentum conservation. 
Since the stress scales as a ~ "qv/r, we have v ~ 1/r/. This 
explains the logarithmic size dependence of the diffusion 
coefficient. 

The dynamic structure factor is a quantity readily ac- 
cessible through scattering experiments. For small wave 
numbers R g k -C 1, only the center of mass motion of the 
polymer can be captured. When R g k 3> 1 and t is much 
less than the relaxation times, the screening lengths be- 
come important. In the limit of 6 ^> 1, S' SD (k, t) shows a 
stretched exponential decay with an exponent 2/3, and so 
does S' ES (k, t) for e ^ 1 with an exponent 1/2. Again the 
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role of the outer environment is reflected in the decay rates 
r SB and r ES given by eqs. (l54l) and fl70l) . respectively. 
In the free membrane case, the dependence r SB ~ k 3 
resembles that of the Zimm model in 3D. For confined 
membranes on the other hand the behavior _T ES ~ k 4 is 
analogous to that obtained from the Rouse model. 

The dynamics of a hydrophobic polymer embedded in 
a 2D membrane was previously investigated by Muthuku- 
mar both for a Gaussian polymer and a polymer with ex- 
cluded volume effects |20) . In his treatment, the membrane 
was regarded as an isolated entity without any couplings 
to the outer environment. The membrane 2D nature was 
taken into account through an anisotropic viscosity. It was 
shown that the longest relaxation time is proportional to 
p^ 1 for a Gaussian chain. This agrees with our limiting 
expressions of the relaxation times in eqs. (1451) and (|61D 
obtained for small <5 and e, respectively. In ref. [3D], the ex- 
cluded volume effects were taken into account by a mode- 
dependent polymer blob size. In this case, he showed that 
the relaxation time scales as p~ 3 / 2 . As discussed in sect.|5J 
we find the mode dependencies of the relaxation times are 
altered in the presence of excluded volume interactions. In 
the small size limit, we indeed recover the scaling p~ 3 l 2 for 
the free and confined membrane cases (see table [T]) . No- 
tice again that the polymer dynamics in this limit remains 
unaffected by the outer environment. For large polymer 
sizes, we obtain either t sd ~ p~ 9 / 4 or r ES ~ p~~ 3 for real 
polymer chains. 

A related situation to a polymer in a confined mem- 
brane is a dilute polymer solution trapped in a slit-like 
geometry whose width is much smaller than the polymer 
blob size. Using scaling arguments, Brochard calculated 
the polymer relaxation time scales as p~ 5 / 2 [31]. An ex- 
perimental realization of such a geometry was done by Lin 
et al. who confined dilute DNA solution in quasi-2D 110 
nm wide slits |34) . The relaxation times measured in this 
case was found to scale as p~ 22 . However, it should be 
emphasized that this scenario is different from the model 
discussed in this article. In our model, the polymer chain 
is strictly confined to the 2D plane of the membrane which 
itself is embedded in a 3D bulk fluid. 

At this stage, a rough estimate of the screening lengths 
would be useful. As reported in ref. [2j , the membrane vis- 
cosity of dimyristoylphosphatidylcholine bilayers at 32°C 
(rounded to the nearest order) is approximately 0.1 Ns/m 2 
and the viscosity of water is rj s w 10~ 3 Ns/m 2 . For sup- 
ported membranes we can approximate the height of the 
intervening solvent region to be h « 10~ 8 m [Mj. Hence 
we obtain v^ 1 w 2.5 x 10~ 7 m and « 0.5 x 10~ 7 m. As 
described in the Introduction, the experiments with DNA 
adsorbed on supported membranes point towards a Rouse- 
like behavior [2l)ll27| . The strong electrostatic attraction 
between the negatively charged DNA and the positively 
charged membrane prevents any out-of-plane motions of 
the polymer. Because the typical length scale of a DNA 
molecule used in the experiments were of several microns, 
the scenario should be close to the case of e s> 1. Hence 
our result is consistent with the experimental observations 
showing the Rouse-like behavior. A more recent experi- 



ment on DNA adsorption on free standing cationic giant 
unilamellar vesicles showed that the diffusion coefficient 
of DNA molecules lies in the crossover region between the 
logarithmic to the algebraic regimes [29] . 

We finally discuss the mobility tensor for a supported 
membrane. In this membrane sits at the bottom 

of a trough filled with a bulk solution having an infi- 
nite depth when compared to the membrane thickness. In 
fig. [U this corresponds to the case where h + is infinitely 
large while h~ is finite. Using the general expression of 
the mobility tensor in eq. ©, we obtain 



G a p [k] 



1 



rjk 2 + k[r]s + T] s coth(kh~)] 



5 at 



k a k 



k 2 



(80) 

If we further assume that h is very small, the above 
equation reduces to 



G Q/3 [k] = 



1 



rjk 2 +r]tk + (rj s /h ) 



Sal: 



k a kf 



(81) 



This mobility tensor for a supported membrane contains 
two screening length scales, i.e., 77/77+ and (rjh~ /r/^) 1 ^ 2 . 
The investigation using the above mobility tensor is left 
as our future work. 
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Appendix A. Derivation of the general mobility 
tensor 

Our derivation of the mobility tensor for the membrane 
closely follows that by Inaura and Fujitani [35 . As pre- 
sented in fig. [2j we consider a more general case where 
the two walls are located at different distances from the 
membrane, i.e., h + ^ h~ . The membrane is assumed to 
be impermeable. 

Our purpose is to derive the in-plane force f s on the 
membrane due to the bulk solvent and walls (see eq. @). 
We first take the Fourier transform of v ± (r, z) by 



dVv ± 



(r, z) exp(— ik ■ r), 



(A.l) 



where r = (x,y) and k = (k Xl k y ). The projection of 
the vector v+ [k, z] on the xy-p\&ne can be expressed as 
[k, z]k + tr^fk, z]k where k = (k x /k,k v /k) and k = 



(—k y /k,k x /k) with k — |k|. The subscripts || and _L indi- 
cate the components parallel and perpendicular to k, re- 
spectively. From eq. (j4|) , the vertical component [k, z] 
obeys the equation 



(A.2) 
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Fig. 6. Projection of v^k, z] onto the icy-plane which is de- 
composed as wj^k + trfk. 



The solution to the above equation can be written as 



■ ± [k,z]=A 1 e- kz +A 2 e kz , 



(A.3) 



with unknown coefficients A\ and A 2 . These coefhcients 
are determined by the stick boundary conditions imposed 
at the membrane-solvent and solvent- wall boundaries: 



where Uj_[k] = v[k] • k and see eq.((7]) for the definition of 
v[k]. Thus we have 

^[k,z] = T^[k] Sinh(fc ^,T + ^ )) . (A.5) 



sinh(fc/i ± ) 



From eq. (|A.5|i . the in-plane force on the membrane 
due to the outer solvent and walls becomes 

/ sX [k] = f s [k] • k 

Q 

= —rj^v±kcoth(kh + ) — i]~v±kcoth(kh~). (A. 6) 

Since the incompressibility condition of the membrane 
fluid implies i>ii[k] = v[k] • k = 0, the perpendicular com- 
ponent of the Fourier transform of eq. @ gives 



7?A;V[k] + / s x[k]+Fx[k] = 0, 



(A.7) 



with -Fj_[k] = F[k] • k. Hence the mobility tensor dehned 
by v[k] = G[k] • P[k] is given by eq. ©. 



Appendix B. Mobility tensors in real space 
Free membrane case 



The mobility tensor in Fourier space is given by eq. (|11[) . 
The expression for G^2 (r) can be found by assuming 



G^(r) = B 1 5 af} +B 2 



(B.l) 



with two coefficients B\ and B 2 . By considering the di- 



agonal and off-diagonal parts of eq. (jB.ll) separately, we 
have 



2B 1 + B 2 = 



1 

2ixr\ 



dfc 



Jo(kr) 



and 

B\ + B 2 





1 




4r? 


1 

27T7? 




/ dfc- 

Jo 


1 


POO 

/ dfc- 

Jo 


1 


2 


4/7 


TTV 2 r 2 



k + v 
Ji(kr) 



(B.2) 
(B.3) 

(B.4) 
(B.5) 
(B.6) 



See ref. [55] for the integral of eq. (IB. 3ft . In evaluating the 
integral of eq. (IB. 61) . we have made use of the following 
relations; 



kr{k + v) 
Jo(fcr) + J 2 (kr) 
k + v 
, H x (i/r) Yxiyr) 



.h(z) 



d 2 



n dzJo(z) d^7T^ = d^J 



dz 



z + a 

Solving eqs. (|B.3[) and (IB. 61) for Si and B 2 , we get 

2 



(B.7) 

(B.8) 



Bi = — 



1 

4r? 



Hi(it) , Fi(^r) 



So 



1 

4r; 



2Hi(i/r) 



2Yi(vr) 



-n (ur)+Y (vr) 



(B.9) 



(B.10) 



The preaveraging of eq. (|B.1|) yields 

5 SD (n-m) = i {2B X +B 2 ) 
Hence we obtain eq. (173")) using eq. (IB. 3ft . 



(B.ll) 
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Confined membrane case 

Now the mobility tensor in the Fourier space is eq. (fT4|) . 
Similar to the free membrane case, its real space expres- 
sion can be written as 



(B.12) 



with two coefficients C\ and Ca. Then it follows that 
1 



and 



Ci +C 2 



27TT7 Jo r(fc 2 + k J 

1 I" 1 JCi(/ir) 
27r?7 K 2 r 2 Kr 



Solving these equations, we obtain 



1 1ixr\ 



K (Kr) + 



Ki(kt) 



1 



C 2 = 



-K { K r) - 



2K 1 (nr) 



(B.13) 
(B.14) 

(B.15) 
(B.16) 

(B.17) 
(B.18) 
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